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Abstract 

The first measurement of tlie radiative muon capture (RMC) rate on a proton was 
recently carried out at TRIUMF. The TRIUMF group analyzed the RMC rate, 
r^Q, in terms of the theoretical formula of Beder and Fearing, and found the sur- 
prising result that gp = fp{(f' = — 0.88m^) is 1.5 times the value expected from 
PCAC. To assess the reliability of the theoretical framework used by the TRIUMF 
group to relate Frmc to the pseudoscalar form factor fp, we calculate Frmc in 
chiral perturbation theory, which provides a systematic framework to describe all 
the vertices involved in RMC, fulfilling gauge-invariance and chiral-symmetry re- 
quirements in a transparent manner. As a first step we present a chiral perturbation 
calculation at tree level which includes sub- leading order terms. 

Key words: fj,~ — > 1/71,7, heavy baryon chiral perturbation theory, pseudoscalar 
coupling 
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1 Introduction 



It has long been a great experimental challenge to observe radiative muon 
capture (RMC) on the proton, /i~ + p — > n + z/^ + 7, because of its extremely 
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small branching ratio. Recently, an experimental group at TRIUMF [1] finally 
succeeded in measuring Frmc; the capture rate for RMC on a proton|^ The 
matrix element of the hadronic charged weak current h'^ = — between 
a proton and a neutron is given by 



n{pf)\V'-A'\p{p,^) 



u{pf) 



fv{q n + -7^ ^ + fA{q n 75 + q 75 



u{pi 



where q = Pi — Pj, and the absence of the second-class current is assumed. 
Of the four form factors appearing in eq.(l), fp is experimentally the least 
well known. Although ordinary muon capture (OMC) on a proton, /i~ + p — 
n + P/j,, can in principle give information on fp, its sensitivity to fp is intrin- 
sically suppressed. This is because the momentum transfer involved in OMC, 
= — 0.88m^, is far away from the pion-pole position = m^, where the 
contribution of fp{q^) becomes most important. RMC on a proton provides a 
more sensitive probe of fp than OMC, because the three-body final state in 
RMC allows one to come closer to the pion pole. 

To relate Frmc to /p, the authors of [1] used the theoretical framework of 
Beder and Fearing [2]. In this framework, as in many earlier works [3-6], one 
invokes a minimal substitution to generate the RMC transition amplitude from 
the transition amplitude for OMC, the hadronic part of which is given by Eq. 
(1). The actual procedure used in [2] is as follows. First, the pion-pole factor is 
exphcitly extracted from fp as fp{q^) = fp/ {q"^ — nT-l), where fp is a constant. 
Then one replaces every q in Eq.(l) with q—eA {A is the electromagnetic field) 
except the q appearing in the q"^ dependence of fv, fA and /m- ^kmg resulting 
from this treatment has a parametric dependence on fp. In the analysis in 
ref. [1], fp is adjusted to optimize agreement between F^^^^ and the measured 
rate F^j^^ (more precisely -R(> 60 MeV)). The result of this optimization, 
expressed in terms of gp = fp{q^ = — 0.88m^) = /p/(— 0.88m^ — m^), is 
gp = (10.0 ± 0.9 ± 0.3)gA, where gA = /yi(0). This value is ~ 1.5 times the 
value expected from PCAC. This surprising result should be contrasted with 
the fact that gp measured in OMC is consistent with the PCAC prediction 
within large experimental uncertainties [7]. 

A natural question one could ask is: How reliable is T^^q used in deduc- 
ing gp from F^jJJq ? It seems important to reexamine the reliability of the 
existing phenomenological approach [2] which uses a selective minimal substi- 
tution. Chiral perturbation theory (ChPT) provides a systematic framework 
to describe the electromagnetic-, weak-, and strong-interaction vertices in a 



To be more precise, the TRIUMF experiment determined the partial capture 
rate R{> 60MeV), corresponding to emission of a photon with > 60MeV. 
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consistent manner, thereby allowing us to avoid applying a phenomenolog- 
ical minimal-coupling substitution at the level of the transition amplitude. 
Furthermore, ChPT enables us to satisfy the gauge-invariance and chiral- 
symmetry requirements in a transparent way. 

Starting with the seminal work of Gasser and Leutwyler [8] ChPT has proven 
to be a very powerful and successful technique for hadronic phenomenology 
at low energies [9]- [12]. Muon capture is another favorable case for applying 
ChPT since momentum transfers involved here do not exceed m^, and is 
small compared to the chiral scale A ~ 1 GeV, indicating the possibihty of 
a reasonably rapid convergence of the chiral expansion. In the case of OMC, 
Bernard et al.[14] and Fearing et al.[15] used heavy-baryon ChPT to evaluate 
fp with better accuracy than achieved in the PCAC approach. In the case of 
RMC, a ChPT calculation provides a natural extension of the classic work of 
Adler and Dothan[13] based on the low-energy theorems. These observations 
motivate us to attempt a systematic ChPT calculation of Prmc- As a first 
step we calculate the total capture rate Frmc and the spectrum of the emitted 
photons, dr^Mc{k)/dk, to sub-leading order in chiral perturbation expansion. 
Thus, our calculation includes nucleon recoil contributions of 0{1/M). We 
hmit ourselves here to the case of RMC from the fi-p atom with statistical spin 
distributions, leaving out the hyperfine-state decomposition and the treatment 
of RMC from the p^p molecule. 



2 Calculational Method 

We employ heavy-baryon chiral perturbation theory [16] and use the effective 
Lagrangian C^h as given in [10]. Cch is written in the most general form in- 
volving pions and heavy nucleons in external weak- and electromagnetic-fields 
consistent with chiral symmetry. We expand Cch in increasing chiral order as: 

£eh=4°^+£S + 4'i + --- ■ (2) 

Here C^"^ represents terms of chiral order u given by z/ = c? + |n — 2, where d 
is the summed power of the derivative and the pion mass, and n denotes the 
number of nucleon fields involved in a given term [17]. We limit ourselves here 
to a next-to-leading chiral order (NLO) calculation and therefore we only keep 
terms with z/ = and z/ = 1. To this chiral order we need only consider tree 
diagrams, and then C^^lf simply represents 1/M "nucleon recoil" corrections 
to the leading "static" part C^lf. We give below the explicit expressions for 
the Cflf and C^^lf, in which only terms of direct relevance for our NLO 
calculation are retained. 
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jrio) ^ J^^^ [D^UD^U] + ... 
= N{^v■D + gAS■u}N 



(1) 

■kN 



I 9 A, 



2M 



(v-D) 



D. D- -^{S -D^vu}^ 

2M 2K'r ' ^ 



AM 



[S^, ((1 + K.)Tr/+ ) }iV + . . . 



(3) 
(4) 

(5) 



Here U = yl — tF^//^ + ■ vr//^ denotes the chiral field in the sigma gauge, 
and the heavy nucleon spinor of mass M. We have also used other standard 
notations, see [10]: 



F^^v, + A„ f;:^v,-a, 

f;.^u^F^,u + uF;:,u^ . (6) 

The covariant derivatives above include the external vector and axial vector 
fields, V/_j = and A^ = ^^-y, respectively. If we choose the four- velocity 
to be f'^ = (1,0), the spin operator of the heavy nucleon becomes 
= (0, |(t). The only parameters appearing in the above expressions are the 
pion decay constant, = 93 MeV, the axial vector coupling, = 1.26, and 
the nucleon isoscalar and isovector anomalous magnetic moments. Kg = —0.12 
and = 3.71. Thus, to the chiral order of our interest, Cch is well determined. 




Fig. 1. 



We consider all possible Feynman diagrams up to chiral order i/ = 1 which 
contribute to the process fj,~+p^n + h' + 'j. These are displayed in Figs.l- 
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Fig. 2. 




Fig. 3. 




Fig. 4. 



6. The zigzag lines in these diagrams represent the W~ boson that cou- 
ples to the leptonic and hadronic currents in the standard manner. In the 
actual calculation, taking the limit mw oo, we make the substitution: 
W~ — >• (V~ — A~){t^ — ir^)/2, and treat V and A as static external vector 
and axial sources, respectively. Then the diagrams in Figs. 1-6 reduce to those 
that would result from the simple current-current interaction of the V — A 
form. The reason for explicitly retaining the W~ boson lines is to clearly sep- 
arate the different photon vertices (see e.g. Fig.6). The leptonic vertices in 
these Feynman diagrams are of course well known. The hadronic vertices are 
obtained by expanding the ChPT Lagrangian [Eqs. (2), (3), (4) and (5)] in 
terms of the elementary fields N, n, V and A and their derivatives. The leading 
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Fig. 5. 




Fig. 6. 



order terms arise from whereas the NLO contributions are 1/M "recoil" 
corrections due to d^^lj- The evaluation of the transition amplitudes corre- 
sponding to these Feynman diagrams is straightforward. We denote by Mj 
(i = 1 . . . 6) the invariant transition amplitudes corresponding to Fig.(l)-(6), 
respectively. They are given by: 



Ml = e^(A) 



Mi 



2{k ■ /i) 



Hi{a)hlH,{a')] (7) 



2,3,4,5,6, (8) 



which include the following hadronic operators: 



hi 



{qiY - ml 



+ ^v'^ S ■ ij) + n) — ■ {p + n) 



(9) 
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2M Ep - Uk 



[v^ - 2gAS^) + 2gA 



(qn)^ - ml 



e(A) • (2p -k) + {2 + Ks + K^){-i)e^p^pe'^{X)k^v''SP 



(10) 



hl{X) 



1 



2M + cuk 



{v^ - 2gAS^) + 2gA 



{qNy - ml 



(11) 



hl{X) 



{qMY{2qL + k)-t{X) 
{q% - ml){ql ~ ml) 

'^9a{S ■ qL) - ■ {p + n){qL ■ v) 



[12) 



hliX) = 2 gA 



(gjv)^ - ml 



{S-e{X)) + 



9A {q 



•n 



M (g^)2 - ml 



[v qM){S ■ e{X)) 



M^^ ■ '^^^^""^ " ^ + 2F^^ + K^)ie-^''^^eMvpS, (13) 



^6(A) 



ql-ml 



2gA{S ■ qL) - ■{p + n){v- qL) 



(14) 



In these expressions, ix^u^p = {Ep,p), n = {En, n) and k — {cuk, k) are the four- 
momenta of the muon, neutrino, proton, neutron and photon, respectively. 

The 2;-componcnts of the spins of the muon, neutrino, proton and neutron 
are denoted by s, s' , o' and cr, respectively, while e(A) stands for the photon 
polarization vector. We have also defined qL = n — p and q^ = n — p + k. 

The pion-polc diagrams. Figs. 1(b), 2(b), 3(b), 4, 5(b) and 6, originate from 
£^°\ Eq.(3). The coupling of the axial vector to the tt generates these Feynman 
diagrams. In ChPT the pion-pole contributions, which arise automatically 
from a well-defined chiral Lagrangian, are completely determined by the chiral 
Lagrangian. The fact that they need not be put in by hand constitutes a 
major advantage of the ChPT approach over the phenomenological approaches 
which have been used in the earlier calculations [2-4]. For example, the term 
originating from Fig. 5(b) does not appear in Ref. [4]. In addition, due to the 
pure pseudoscalar pion nucleon coupling, the pion-pole terms are proportional 
to 1/M in Ref. [4]. In this context it is also worthwhile to mention that the 
pseudoscalar coupling gp itself does not appear exphcitly in ChPT calculations 
of the transition amplitudes since gp is effectively accounted for via the pion- 
pole diagrams. As mentioned in the introduction, Cch determines gp [14,15]. 
However, since the same Cch directly determines the transition amplitude of 
RMC, gp does not feature in our expressions for Mj's. 
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It is safe to assume both the muon and the proton to be at rest by neglecting 
the binding and kinetic energies of the fip atom. Thus, /i = [m^, 0) and p = 
(M, 0). For the neutron four-momentum n, we retain its three- momentum 
n but neglect the recoil energy, oi = M + fP/2M ^ M. The maximal 
value of \n\ equals m^^ giving a recoil energy n^/2M 6 MeV, which is small 
even compared with m^. With n ~ {M,n), we have ql — (0, n) and qn — 
{ujk,n + k). Consequently, all terms proportional to v -qi vanish. We choose to 
work in the Coulomb gauge with the result f ■e(A) = 0. With this gauge choice 
and the above kinematical approximations, the hadronic radiation diagrams. 
Figs. (2) and (3), become 0{1/M^) [see Eqs. (10) and (11)], and therefore do 
not contribute to the chiral order under consideration. Moreover, in the sum 
^2 + ^3) the terms proportional to vanish in our approach (to the order 
under consideration), whereas in the treatment of, e.g., [3], these terms are 
numerically large. 



3 Numerical Results 



As stated, wc consider here only the RMC from the /xp atomic state with 
the hyperfine states unseparated. Within our kinematical approximations the 
spin-averaged total capture rate is given by 



trmc = (^^J \m\' H^^r J 



cPn r d^u r d^k 1 



(27r)3 J (27r)3 J (27r)3 2^;^ 



aa'ss'X 



where the sum is over all spin and polarization orientations, M = X]i=i 
with Mi given by Eqs. (7) and (8); $(0) is the value of the fip atomic wave- 
function at the origin. In the kinematical approximation stated earlier, Eq.(15) 
simplifies as 



('^fc)max~»Tl^ 

Trmc = (87r^C) j duJkUJk{m^- uJkf 



X [dcose y: E (m.m;) , (16) 

ij=l,4,5,6 cra'ss'X \^\=mij,-u,k 

where we have introduced the abbreviation C = (eC/y^)^ (1/2^71^). The eval- 
uation of the spin sum is tedious but straightforward; the resulting lengthy 
expressions will be given elsewhere. 
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Table 1 summarizes the mimcrical values for the total capture rate Frmc- 
We show in the table the breakdown of Frmc into the leading- order contribu- 
tion (9((1/M)°), coming from C^^lf, and the next-to-leading-order contribution 
0{1/M), arising from We also show the value of Frmc which would result 
if all the diagrams containing the pion-pole, Fig. 1(b), 4, 5(b) and 6, are omit- 
ted. Our result for the total capture rate Frmc = 0.075s~^ is close to the value 
given in [4], Frmc = 0.069s~^, and practically identical to Frmc = 0.076s~^ 
reported in [5]. Our 0{1/M) recoil corrections account for about 20% of the 
leading order 0((1/M)°) contribution, which indicates a reasonable conver- 
gence of the chiral expansion. It should be noted that the size of the 1/M 
corrections is noticeably larger in the approach of [5] . As one can see from Ta- 
ble 1, about 30% of the total value of Frmc comes from the pion-pole exchange 
diagrams. 

Table 1 



Total RMC capture rate in s 





Frmc 


Frmc without TT 


O((l/M)0) 


0.061 


0.043 


0{1/M) 


0.014 


0.010 


total 


0.075 


0.053 



In Fig. 7 we plot the spectrum of the emitted photons dT RMci^-^ k) / duik- In ad- 
dition to the result of the full calculation, the figure includes the spectrum cor- 
responding to the leading-order calculation, i.e., the 0{{1/M)^) contribution 
only. For the sake of comparison, we also show the result of [2,5] corresponding 
to the use of the Goldberger-Treiman value gp = 6.651^. 



4 Discussion and Conclusions 

A direct comparison of our calculation with the experimental data [1] is pre- 
mature because we have not considered capture from the singlet and triplet 
hyperfine states separately, or capture from the p/ip molecular state. This also 
means that at this stage we cannot directly address the "gfp problem" that 
arose from the TRIUMF data [1]. However, it is worthwhile to make the fol- 
lowing remark. As one can see from Fig. 7 our ChPT calculation gives for the 
spin-averaged /^.p-atomic RMC a photon spectrum that is slightly harder (by 
about 10% for > 60 MeV) than what was obtained in [2,5] with the use 
of the G-T value, gp — 6.651^. Meanwhile, as mentioned earlier, ChPT gives 
a value of gp consistent with gp = 6.6gA [14,15]. Thus, there is a possibility 
that, even with the same value of gp, a ChPT calculation gives a somewhat 
harder 7 spectrum than the conventional method. It remains to be seen to 
what extent such a difference in the 7 spectrum influences the gp value de- 
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duced from the experimental spectrum in the higher energy region. Of course, 
a more quantitative statement can be made only after a more detailed ChPT 
calculation becomes available in which the hyperfine states are separated and 
the p/xp-molecular absorption is evaluated. We also remark that using relativis- 
tic kinematics, instead of the kinematic approximation employed in Eq.(16), 
softens the photon spectrum to a certain extent. 

We repeat that the present calculation includes only up to the next-to-leading 
chiral order (NLO) contributions. The next-to-next-to-leading order (NNLO) 
calculations are obviously desirable. For this one must include the V — 2 chiral 
Lagrangian, and C^^l^, and also loop corrections arising from C'^^ and 
>C^°^. The finite contributions from the loop diagrams would give momentum- 
dependent vertices, which would correspond to the form factors in the language 
of the phenomenological approach [2,3,5]. These contributions are probably 
small but it would be reassuring to check that explicitly. One problem in 
extending the present calculation to the next order is that, although the forms 
of vC^^-* and d'^li have been determined [8,18], some coefficients of the counter 
terms in C\l^ still remain undetermined. On the other hand, chiral expansion 
for muon capture is characterized by the expansion parameter m^/M, and is 




Fig. 7. Spectrum of the emitted photons. The full line represents the full calculation 
including 0{{1/M)^) and 0((1/M)^); the dashed line represents the result that 
contains only the 0{{1/M)^) contributions; the dotted line shows the result of [2,5] 
with gp = Q.&gA- 
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expected to converge reasonably rapidly. Indeed, in the case of OMC, where 
the V = 2 calculation is much less involved, explicit evaluations [14,15] show 
that the NNLO contributions amount only to a few percents. It is likely that, 
in the case of RMC as well, NNLO corrections modify our results only by a 
few percents. In this connection we also note that the formalism of Bernard et 
al.[10] used here does not contain the explicit A degree of freedom in contrast 
to the approaches of [16]. Although it is desirable to examine the importance of 
the A, we relegate that to future studies. [After the completion of the present 
work we learned of the first attempt at an NNLO calculation by Ando and 
Min [19].] 
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